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I. INTRODUCTION 

In chaotic and. weakly disordered quantum systems 
which are not self-averaging, phase coherence gives rise 
to sample-to-sample fluctuations in most transport prop- 
erties with respect to a small perturbation in the incident 
energy, an applied magnetic field or the shape of the sys- 
tem. Those fluctuations are universal pi 0] and depend 
only on the symmetry properties, such as the presence 
or absence of time reversal invariance (TRI), and spa- 
tial symmetry 0, 0, 01, El An statistical analysis is well 
described by random matrix theory (RMT) Q • 

The parametric dependence of the conductance has 
been studied experimentally by considering ballistic 
quantum dots connected to electron reservoirs by bal- 
listic points contacts with few propagating modes 0, 0, 
IH IIIllIl- RMT predictions can also be verified in wave 
scatterin g ex perimental systems, such as microwave cav- 
ities [l3l acoustic resonators 0, or elastic media 
|l6j |. where the external parameters are easy to control. 
However, absorption is always present in these experi- 
ments and its influence on the universal transport prop- 
erties is rather dramatic [l^; therefore, many theoretical 
and experimental works have been devoted to the effect 
of absorption on the transmission T and reflection R co- 
efficients of the cavity H E |H EJ |H i| . The 
derivative of those coefficients with respect to the exter- 
nal parameter has not been considered in the presence of 
absorption. A parametric derivative is very important in 
the characterization of mesoscopic systems with a chaotic 
classical limit uA |25| . since it is analogous to the level 

velocity in EalHlia . 

Motivated by recent experiments in microwave cavi- 
ties 0, |23| , in the present paper we study the statisti- 
cal fluctuations of the parametric derivative of T and R 
with respect to an external parameter X, dT/dX and 
dR/dX, in the presence of absorption. We consider a 
chaotic cavity connected to two waveguides with an arbi- 
trary number of channels, with and without TRI, and we 
address both asymmetric and left-right (LR) symmetric 
cavities. As an external parameter we will take shape de- 



formations. The purpose of this work is three fold: first, 
the calculations here presented help to understand the 
distribution of the energy derivative of T in the presence 
of absorption; in fact, we now present a complete theo- 
retical derivation of some of the results used in Ref. |3(]. 
Second, they also can serve to motivate the experimen- 
tal analysis of the distribution of the derivative of T but 
with respect to shape deformations, where the results of 
the present paper can be applied. That is the case of 
Ref. |33 where, in order to improve statistics, the shape 
is modified by varying one lenght of the resonator used 
in the experiments. Fin ally, in a similar way, the experi- 
mental situation of Ref. |23| can be used as well to study 
energy and shape deformation derivatives of R. 

The results presented here are valid for strong absorp- 
tion. However, they reproduce those existing in the lit- 
erature for the distribution of dT/dX at zero absorption 
intensity p3. l25j . In the absence of absorption the dis- 
tribution of the parametric conductance derivative was 
calculated analytically by Brouwer et al |24[ for an asym- 
metric quantum dot with two single- mode point contacts. 
The 9T/<9AT-distribution has algebraic tails and in the 
absence (presence) of TRI it shows a cusp (divergence) 
at zero derivative; the second moment is finite (infinite). 
The reflection symmetric case was considered in Ref. l2a 
There, the distribution of dT/dX diverges logarithmi- 
cally at zero derivative, it has algebraic tails with an ex- 
ponent which is different to that of the asymmetric case. 

The paper is organized as follows. In Sect. [H] we 
present the main formal elements used throughout the 
paper, such as the scattering matrix S and its parametric 
derivative in the presence of absorption. Sect. Ill Al is ded- 
icated to asymmetric cavities. The Poisson kernel for S 
and its application to chaotic scattering in the presence of 
absorption is presented by means of a phenomenological 
model; the parametric derivative of S is defined in terms 
of a Wigner time-delay matrix whose eigenvalues are the 
proper time-delays, the inverse of them being distributed 
according to the Laguerre ensemble. The general struc- 
ture for S and its parametric derivative for cavities with 
LR symmetry is introduced in Sect. IIIBI The mean 
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FIG. 1: A ballistic chaotic cavity connected to two leads with 
Ni, N2 channels. N p equivalent "parasitic" channels are at- 
tached to the cavity by tunnel barriers with transmission T v 
pj^| . The absorption strength is given 7 P = N P T P in the limit 
N p — ► 00, T p — ► keeping the product constant (2^1 . 



and variance of the parametric velocities for T and R, as 
well as the correlator between the channel-channel trans- 
mission and reflection coefficients are calculated in Sect. 
If II Al in the presence of TRI, and in Sect. IIII Bl in its 
absence. Sect. IIVI is dedicated to LR-symmetric cavities 
where we calculate the variances of parametric velocities 
in the presence (absence) of TRI. Finally, a summary of 
the results as well as the conclusions are presented in 
Sect. 



II. THE S MATRIX AND ITS PARAMETRIC 
DERIVATIVE 

A. Chaotic scattering by asymmetric cavities in 
the presence of absorption 

The scattering problem of a ballistic cavity connected 
to two waveguides, each supporting Nt, N2 transverse 
propagating modes (see Fig. [TJ, can be described by the 
scattering matrix S which, in the stationary case, relates 
the outgoing to the incoming wave amplitudes [3l|. 

The absorption in the cavity is modeled attaching N p 
equivalent non transmitting or "parasitic" channels to 
the cavity by means of a tunnel barrier with transmission 
T p for each one p3. l3^ | . The S matrix is N dimensional 
(TV = iVi + 7V2 + N p ) with a structure given by 



Sll Sl2 Sip 
S21 S 2 2 s 2p 

Spl Sp2 Spp 




(1) 



where the set of indices {1}, {2}, {p} label the Ni, N2, 
N p channels. Here, the submatrix S of dimension N1+N2 
describes the scattering problem of the absorbing system. 
The absorption can be quantified by the parameter j p — 
NpTp in the limit N„ — > 00, T p — ► while keeping the 
product constant |24j . 

T and R are obtained from S, actually only from S, as 
follows 



T = \ s *»\ 2 and R = E \ Sa » 



(2) 



In our case only two of the three basic symmetry classes 
in the Dyson's scheme 33] are relevant . We assume that 
S satisfy flux conservation by the restriction 



N ■ 



(3) 



where Ijv stands for the unit matrix of dimension N. 
This case is called "unitary" and it is designated as (3 = 2. 
In addition, in the presence of time reversal invariance S 
is symmetric, 



S = S J 



(4) 



This is the "orthogonal" case, designated as (3 = 1. Note 
that the N p channels are normal scattering channels for 
the matrix S, while they are absorbing channels for the 
matrix S, which is a subunitary one and describes the 
physical system; it represents the scattering matrix of 
the absorbing system where the flux is not conserved. 

For systems with a chaotic classical limit, most trans- 
port properties are sample specific and a statistical anal- 
ysis of the quantum-mechanical problem is needed. That 
study is performed by the construction of ensembles of 
physical systems, described mathematically by ensembles 
of S matrices distributed according to a probability law. 
The starting point is a uniform distribution where S is a 
member of one of the circular ensembles: circular unitary 
(orthogonal) ensemble, CUE (COE), for (3 = 2 {(3 = 1) 

In the presence of direct processes, the information- 
theoretic approach of Refs. |3f|, |3^ leads to an S matrix 
distributed according to Poisson's kernel [33 



Pg\S) = c 



[det(l w -(5)(5)t)] ( ^ +2 -« /2 



|det(ljv-5(5)t)| 



f3N+2-p 



(5) 



where (S) is the ensemble averaged S matrix. 

A useful model to construct the Poisson ensemble con- 
sist of a cavity connected to leads by tunnel barriers [3^| . 
In the case we are concerned with, where only the ficti- 
tious waveguide contains a tunnel barrier, the averaged 
S matrix can be written as 



(S) 




V^Tpt 



(6) 



A',. 



a,b£l 



As before, l n stands for the unit matrix of dimensions n 
and 0„ for the n-dimensional null matrix. 

In what follows we restrict ourselves to the case where 
T p = 1, i.e. Pj!p (S) is just a constant and the S matrix 
is uniformly distributed. In this case, we are restricted 
to a strong absorption situation, where the parameter 7 P 
takes only integer values (j p = N p ). Also, the results 
here presented are valid for no absorption (N p = 0), and 
a simple extrapolation to non integer values of 7 P is qual- 
itatively correct, as will show later on. 

If the coupling to the fictitious waveguide is perfect, 
we can use the well known definition of the parametric 
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derivative of S. The derivative of S with respect to the 
energy of incidence E can be defined in terms of a sym- 
metrized form of the Wigner-Smith time delay matrix 
[3^ | , whose eigenvalues are identical among them [2(j . In 
dimensionless units we have 



dS_ 



1S 1 / 2 Q e S 1 ' 2 , 



(7) 



where we have defined e = 2-kE/ A with A the mean level 
spacing, Q £ is an N X N Hermitian matrix for j3 = 2, 
real symmetric for (3=1. The eigenvalues of Q £ are r^ 1 
times the proper delay times, where tjj = 2irh/A is the 
Heisenberg time. In an analogous way, the derivative of 
S with respect to an external parameter X is defined as 

m 



^-=iS 1 / 2 Q x S 1 '\ 
ox 



(8) 



where we have also defined a dimensionless parameter 
x = X/X c with X c a typical scale for X, and Q x is an 
N x N Hermitian matrix, real symmetric in the presence 
of time-reversal symmetry. 

For classically chaotic cavities the joint distribution of 
S, Q e and Q x is given by [iof 



P (S,Q e ,Q x ) oc(detQ e ) 



-2/3AT-3(l-/3/2) 



x exp 



tr 



(9) 



S is independent of Q £ and Q x , and uniformly distributed 
in the space of scattering matrices. Following ^(jj Qx has 
a Gaussian distribution with a width set by Q £ , that can 
be parametrized as follows |40| 



*- 1T if*-\ 



(10) 



where ^ is a N x iV matrix, complex in the unitary case 
and real in the orthogonal one, such that 



(11) 



and H is a N x N Hermitian matrix for (3 = 2, and real 
symmetric for (3 = 1. H has a Gaussian distribution with 
zero mean and a variance 



(H ab H cd ) = U { 1 



4 5 ad 6 bc 



= 2 



idhc + SacSbd) P = 1 



(12) 



as can be seen by substituting (|10f> and (|llfl into JSJ. 
Now, we diagonalize Q e , 



(13) 



The elements {r n } (n = 1, ...,N) of f are the di- 
mensionless delay times. Their reciprocals x n = l/r n 
(n = 1, . . . , N) are distributed according to the Laguerre 
ensemble |40j . 



a<b 



(14) 




FIG. 2: A ballistic chaotic symmetric cavity connected to two 
leads supporting Ni channels. N p non transmitting channels 
are attached to the cavity to model the absorption. 



The matrix of eigenvectors, W, is uniformly distributed 
in the unitary (orthogonal) group for (3 = 2 (/3 = 1). 

For the calculations we are interested here, it is also 
convenient to parametrize the S matrix and its paramet- 
ric derivative as |4j 

S=UV, ^=iUQ E V, ^ = iUQ x V, (15) 
where U, V are the most general N x N unitary matri- 



ces in the unitary case {(3 = 2), while V 
orthogonal one {(3 = 1). 



U 1 



the 



B. Chaotic scattering by symmetric cavities in the 
presence of absorption 

For a system with spatial left-right (LR) symmetry, 
as shown in Fig. |3 the S matrix is block diagonal in a 
basis of definite parity with respect to reflections, with a 
circular ensemble in each block 0, . 

In the presence of absorption the S matrix that de- 
scribes the scattering of LR ballistic cavity connected to 
two waveguides, is of dimension N = 2Ni + N p , where Ni 
are the number of channels in each waveguide (the two 
waveguides have the same number of channels and are 
symmetrically positioned); N p is the number of absorp- 
tion channels that we assume symmetrically distributed 
in the cavity. In this case, the general structure for S is 

s 



S 



t' 



(16) 



where r', t' are N' x N' matrices, with N' = N x + N p /2. 
They represent the reflection and transmission matrices, 
respectively, associated to the total S matrix given by 
(|16[) . and not for the physical one. The iVi x iVi trans- 
mission and reflection matrices, t and r, associated to 
the system with absorption, are submatrices of t' and r', 
respectively. 

iS-matrices of the form given by Eq. (|ltj|) . which also 
satisfy © are appropriate for systems with reflection 
symmetry in the absence of TRI. With the additional 
condition |0J it is appropriate for LR-systems in the pres- 
ence of TRI |42j. However, when TRI is broken by a uni- 
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form magnetic field, the problem of LR-symmetric cavi- 
ties is mapped [5j to the one of asymmetric cavities with 
(3=1 with t' replaced by r' . 

Matrices with the structure H16|) can be brought to a 
block-diagonal form |4| 



r ■ _ tdT I Si 

K ° I S 2 



Ro , 



(17) 



(18) 



where Rq is the rotation matrix 

_ J_ / t N , 1 N , 

l ra denotes the nxn unit matrix; Si = r'+t', S% = r' — t 1 
are the most general N 1 x N 1 scattering matrices. They 
are statistically uncorrelated and uniformly distributed: 
CUE 03 = 2), COE 03 = 1) 0. 

The transmission and reflection coefficients T and R, 
for LR-symmetric ballistic cavity in the presence of ab- 
sorption are then given by 

JVi 

\(Si) ab -(S 2 ) ab \ 2 and (19) 



/ - y • 1 

a, 6=1 
1 ^ 

R = ^ ^ |(»5'l)a6 + (<S2)afc| 2 , 
a, 6=1 



(20) 



respectively. 

The parametric derivative of S is defined through the 
parametric derivatives of Si and S2 as in Eqs. and 
(jHJ. The joint distribution is satisfied for each ma- 
trix Sj (j = 1, 2). Finally, we note that they can be 
parametrized as in Eqs. i|15|) . 

In what follows we calculate the mean and variance of 
dT/dq and dR/dq, where by q we mean e or x. Also, we 
calculate the correlations between the (/-derivative of the 
channel-channel transmission coefficients. 



The ensemble average of dT/dq and dR/dq can be 
calculated if we substitute the parametrization (|15fl into 
Eqs. (|2*2*|) and (J2HJl. In this way, we get expressions in 
terms of twice the real part of products of averages of 
linear expressions in Q q times averages of nonlinear ex- 
pressions in V and /or U (V = U T for [3 = 1). Using 
the results of Ref. [43, the averages with respect to U or 
V are real positive numbers, while ((Q x )ab) — because 
the matrix H of Eq. (|10|) has zero mean; ((Q £ ) a b) is a 
purely imaginary. Then, the results are 



{dT/dq) = = (dR/dq), q = e,x, 



(24) 



as expected because the distributions of dT/dq and 
OR/Oq are symmetric with respect to the zero derivative 

mm. 

The fluctuations require a more sophisticated analysis. 
Let us define the correlation coefficients by 



Q a'b' 



d(Jab da a ' b ' 



dq dq 

The variances of dT/dq and dR/dq are then given by 



(25) 



(dT/dq)' 



(dR/dqf) = 



with 



1 a'b' 



2Re 



a,a'el b,b'E2 

E E 

a,a'£l 6,b'el 



Q Q* dS* b dSg'b' 

^>ab^> a 'b'~^ ^ 



(26) 
(27) 



r,* r,* dS a b dS a 'b' 

ab a ' b '~d~q dq~ 



(28) 



III. MEAN AND VARIANCE OF dT/dq AND 
dR/dq (q = e, x) FOR ASYMMETRIC CAVITIES 

In this section we first calculate the mean of the q- 
derivative (q = e, x) of T and R. Second, we calcu- 
late correlation coefficient between the g-derivative of two 
channel-channel transmission coefficients, from where, fi- 
nally, we can obtain the variance of dT/dq and dR/dq. 
The present section is devoted to asymmetric cavities for 
both (3=1 and (3 = 2 symmetries. 

By convenience we define the probability to go from 
channel b to channel a as 



Cab = \S a b\ 

from Eqs. J2J) we can write 



dT \ ^ \ ^ da a b , 

-Q- q - LL^- and 

H ael be2 y 

dR \ - do ab 
dq ~ ^ 



a,bel 



dq 



(21) 

(22) 
(23) 



where we have written explicitly the elements of the S 
matrix. 

Because of the complexity of the calculations, in the 
rest of this section we will consider the two symmetries 
(3=1 and (3 = 2 in a separate way. 



A. The orthogonal case 



( 1) 

1. The correlator C\ a , b , 



In the orthogonal case, the substitution of the 
parametrization given by Eq. (|15fl . with V = U T , into 
Eq. l|2*S|l gives the result 

JV N 

C « ] tv = 2Rc E E <(Q«MQ<rW> 

a,fl=l a',/3' = l 

N 

x J2 [J{a,(3,c,c)-J(c,c,a,(3)}, (29) 

c,c' — 1 
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where, in order to simplify the expression, we have de- 
fined the coefficients 

7/ a z\ n <rG,oi,bB ,a' c'b' c 

= {U ai U bS U a , a ,U b ,p,U: a UZ p U:, c ,U* b , c ,) . (30) 

The first (last) two places a, (3 (7, S) of the argument 
of J(a, (3, 7, S), refers to the second and fourth positions 
in the upper (lower) indices of M^" b g a ^, which is 
defined by the second line of Eq. (|30|1 . As we can see 
in App. [S] the rest of the indices of the coefficients M 
are not modified in the construction of Eq. (|29|) . Those 
coefficients M were calculated in Ref. ^3 [ see Eq. (6.3) of 
that reference]; we apply those results to our particular 
case in App. ^ The sums with respect c, c' appearing 
in the second line of Eq. I|29|l give the result 

JV 



^2 [J{a,f3,c,c) - J(c,c,a,(3)] 



c,c' — 1 

= -m 



We substitute Eq. 1(51)1 into Eq. $Z§\i and simplifly to 
obtain a result that depends on n\ and 713 — 112- In App. 
|I|we show that — = Nni [see Eq. <|A11|) ] where rii 
is given by Eq. (|A13ll . Then, we write Eq. (|2"9l as 



(32) 



where 



N 



N 



K( q 1} = N Y,{(Ql) aa )- E ((QiUiQM ■ (33) 

a=l u,f3=l 

is given by Eq. (|33|1 with q replaced by e. K x 
can be written in terms of Q e by direct substitution of 
Eq. i|l(J|) into Eq. Q33[) . The average over the matrix H 
is performed taking into account Eqs. (|12|l and l|llfl for 

(3=1; the result is 

N N 
KP = 4(7V-2)]T<(Q 2 ) qq )+47V £ {{Qe) aa {Qe) PP ) ■ 
a=l a,f3=l 

(34) 

Now, we use the diagonal form of Q s , Eq. l(T3"j l. 
becomes independent of the unitary matrix W, and de- 
pends on two eigenvalues of Q e as 

KP = AN(N - 1) [2 (r 2 ) + N (nr 2 )] , (35) 
K?) = iV(iV-l)[( Tl 2 )_( Tir2 )]. (36) 

The remaining averages of the r variables are performed 
by direct integration using Eq. (|14|l for j3 = 1. (rf) 
diverges for N — 1 , while the next four values of TV give 
the general term 



(rl) 



2N\ 



(N — 2) (N + l)\ ' 



(JV-1)! 
(TlT2> = jN + TfV (3?) 



Then, Eqs. (|35|) and l|36|l are written as 

4iVA e , a £ (iV-2)(iV + l)' 



(38) 



are combined to give the desired 
results for the correlation coefficients, namely 



u i a'b' 
U £ a'b' 



~.n\ab 



(39) 



(N - 2)N 2 (N + 1) 2 {N + 3) 
x {2(l + ^)(l + 5f0 
+ (7V + l)(7V + 2)(<5^ fc b '+^'^') 2 

+2^',( ( , b b ' ( ,°, + S a b 'S a b ,) + 2$6 b a ,6$ 
+ S b a S a b 'S a a , + S b a S b b S a b , + } , 



(40) 



(31) where the dependence on the absorption strength 7 P 



N p is through N = Ni+N 2 + N p . 

From Eqs. 139|) and (|40|l we analyze several cases of 
interest. First, a' = a £ 1, b' = b £ 2, give the vari- 
ances (maximal correlations) of the energy and paramet- 
ric derivatives of the channel-channel transmission coef- 
ficient ddab/dq (q = e, x); those are 



{da ab /dxf) = 4N ((do ab /de) 2 ) 
ff) ,» yA 2(N 2 + N + 2) 



(N -2)N 2 (N + l) 2 (iV + 3) 



(41) 
(42) 



We see that for strong absorption, j p = N p ^> Nl, N2, 
they behave as 

([doab/dx] 2 ) ~ lp \ ({do a b/def) ~ 7 ; 4 - (43) 

Second, when a' = a € 1 and b, b' £ 2, but 6' ^ 6, in the 
limit of strong absorption we obtain 



(44) 



that are smaller compared with the variances given by 
Eqs. I)43|) by a factor of 7" 1 . Finally, when all the in- 
dices are different, in the limit of strong absorption, the 
correlator between the parametric derivatives of two dif- 
ferent single channel transmission coefficients behaves as 



/ d<T ab da ab ' \ 




I da ab 


d(T a b' \ 


\ dx dx 




\ de 


de 1 



dOgb d(J a'b 

dx dx 



dogb do a , v 
de de 



(45) 

which are 7 p 2 times the variances. We conclude that for 
strong absorption, up to the order of {{da a b / dq) 2 ) , the 
correlations between the elements da a b/dq, for a G 1, 
b G 2, are very small. Those quantities enter in the con- 
struction of dT/dq [see Eq. Q22[l] and can be treated as 
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N\N 2 uncorrelated variables with the same distribution. 
This is a relevant simplification when the distribution of 
the parametric derivative of the total transmission coef- 
ficient is desired, assuming the one for each da a b/dq is 
known. That is the case of Ref. 30 where the numerical 
evidence shows an exponential decay for Pi(da a b/de), 
P\(dT/de) being calculated in a very straightforward 
manner. Eqs. (|41|l and l|42[l can be used to obtain the 
decay constant as a function of 7 P [30| . 



2. Statistical fluctuations of dT/dq and dR/dq (q = e, x) 

The second moment of the distribution of dT/dq is 
obtained from Eqs. (|39|) and (|40l) by direct subtitution 
into Eq. I|26|l ; we obtian 



{dT/dxf) = mUer/def 



(46) 



(dT/de) 



NN 2 [(N + l)( 7p + 2) + 2NiN 2 ] 



(N — 2)N 2 (N + 1) 2 (N + 3) 



(47) 



For the particular case N\ = N 2 = 1, Eqs. (|46|l and 
(|4"7|l reduce to Eqs. ljlT|) and (|4^|l . respectively. Also, 
when 7 p = N p = 0, which means no absorption, the 
variance of dT/dq diverges. This is in agreement with 
Ref. |2J where the distribution of dT/dq, was obtained 
in the absence of absorption. The distribution has long 
tails and a divergent second moment. This divergence 
is suppressed in the presence of absorption. We also see 
that the divergence of ((dT/dq) 2 ) disappear when N% or 
N2 is larger than one for any absorption strength. 

In similar way, we substitute Eqs. I|39(l and 14()|l into 
Eq. O to obtain 



(dR/dxY) = 4N((dR/dey} (48) 
,JVi(iVi + l)(N - N!)(N - Ni + 1) 



(dR/de) 



(N — 2)N 2 (N + 1) 2 (N + 3) 



(49) 



for N > 2, while ((dR/dq) 2 ) is infinite for N = 1 as 
N ^ Nx or N Ni - 1. When N = Aq or N = Aq - 1, 
((dR/dq) 2 ) = 0. 

Consider the case Aq = 1 and N 2 = (N = 1 + N p ), 
which is relevant to the experimental data of Ref. |U In 
the absence of absorption 7 P = N p = 0, N = Ni = 1 and 
((dR/dq) 2 ) = as expected (R = 1). For j p = N p < 1, 
((dR/dq) 2 ) is infinite. Again, the divergence is suppresed 
for 7 p = N p > 1. 

We recall that our results are valid in the strong ab- 
sorption regime where by convenience we assumed perfect 
coupling (T p = 1) of the N p absorbing channels. The 
absorption strength 7 P — N p takes only integer values. 
However, a simple extrapolation to T p < 1, which means 
arbitrary 7 p = T p N p , works qualitatively well. In Fig. 
Ek) we compare Eq. (gTJ, for N± — N 2 — 2, with the 



results from numerical simulations |44| for T p — 0.025, 
0.05, 0.075, 0.1, 0.125, and 0.15 with N p = 200 which 
give 7 P = 5, 10, 15, 20, 25, 30. 

In the absence of absorption, i.e. N p — 0, the results 
presented here are strictly valid. In this case R + T = 1 
and the distribution of dR/dq is equal to that of dT/dq. 
In particular their variances are the same: it is easy to 
verify that Eqs. JUJl and (gSJ reduce to Eqs. (gHJ an d 
(|47|l . in complete agreement with the results obtained 
directly from the known distribution of those quantities 
in the absence of absorption j24[. The particular cases 
Ni = 1, N2 = 0, and Aq = N 2 — 1 has been explained 
above. Similar conclusions are valid for the unitary case 
and for f3 = 1, 2 for reflection symmetric case below. 



B. The unitary case 



The correlator C, 



(2)'-' 



The unitary case is simpler than the orthogonal 
one. Following the same procedure, we substitute the 
parametrization (|15|) into Eq. (|28() with the result 



(7(2) ab 

1 a'b' 



2 Re 



N N N 

E E E 

a,/3=l q',/3' = 1 c,c'=1 



{(Qq)pa(Qq)a'P'} 



*m%:%m%;$ - ((Q q u(Q q ) a > ( ,)Mr a %Mf b % 



where we have defined 



M a'b',c'd' Ijji jj, jji* jji* \ 
1V1 ab,cd — \ U ab U cd U a'b' U c'd'/ 



(50) 



with U' a unitary matrix that denotes the unitary ma- 
trices U or V of Eq. 115|) . Those coefficients have been 
l43l 



calculated in Ref. |43J and read 



M 



a b ,c d 
ab,cd 



1 



N 2 - 1 
1 



N 



{°a d c d b °d + °a °c °b °d , 



\ ab 



(51) 



We substitute Eq. I|51|) into C q ' and perform the 
sum over the dummy indices, the result is 



(7(2) ab 

1 a'b' 



l-N(5 a a +6 b b ) + N 2 6 a a '5 b b ' 



N 2 (N 2 - l) 2 



ReX?) 



(52) 

where K q 2 ^ has the same form as Eq. (I33|l but with the 

upper index 1 on the left-hand side replaced by 2, and the 

(2) 

matrix Q q is an Hermitian one. Again, K e is obtained 

(2) 

by replacing q — e. To write K x in terms of Q e we use 
Eq. H10|) and perform the average over H using Eq. Ijl2|l 
for f3 = 2. The result is 

N N 

^ 2) = - 4 E((^)-)+ 4iV E ((^U(QeW). 
a=l a, (3=1 

(53) 



Now, we substitute Eq. I|13fl an d perform the average 
over W to obtain 

= AN(N ~ 1) [(r 2 ) + N ( Tl r 2 )] , (54) 
K?) = N(N-l)[(r 2 )-(r 1 r 2 )] (55) 

By direct integration of the first N terms, Eq. 114fl for 
(3 = 2 give 



2N(N - 2)1 



t it 2 



(N-l)\ 

(N + iy: 



(56) 



(N + l)\ 
Eqs. (52}, O and (5SJ give 

K {2) = , Xi 2) = 1 

Finally, we combine Eqs. 1521) and (|57|l with the result 



(57) 



(-(2) a " 
3S a'b' 



(-(2) ah 
°e a'b' 



l-N(S a a +S b b )+N 2 S-'5 b b 



N 2 (N 2 - l) 2 



(58) 



(59) 



Two different particular cases are of interest. The first 
one, a correlated case, is obtained for a' = a € 1 and 
6' = 6 6 2, for which one obtains that 

(da ab /dxf) = 4N ({da ab /def) , (60) 

\ 2 
(c^/cte) 2 ) = jv2(jV+1)2 , (61) 

which for strong absorption they have the behavior given 
by Eq. (|43[l . Second, uncorrelated cases are obtained 
when a' = a € 1, b' ^ b (6,6' G 2), and when all the 
indices are different, in the strong absorption limit. For 
large j p = N p Eqs. I|44|) and l|45|l are also satisfied for 
(3 = 2. Those quantities are very small compared to 
the order of ((da a b / dq) 2 ) , meaning that in this limit the 
quantities d<j a b/dq for a € 1 and 6 G 2, can be treated 
as -/ViiV 2 uncorrelated variables with the same distribu- 
tion P2{d<J a b/dq). Numerical evidence [3(j also shows an 
exponential decay of -P 2 (<9er a b/cte) for strong absorption; 
the decay constant depends on j p and can be obtained 
from the variance of da a b / dq. 

2. Fluctuations of dT/dq and dR/dq (q — e, x) 

The statistical fluctuations of the energy and paramet- 
ric derivative of the total transmission coefficient is ob- 
tained by direct substitution of Eqs. and l(59"|) into 
Eq. (23) for (3 = 2. The results are 

(dT/dx) 2 ^ = 4N((dT/de) 2 ^ (62) 
2N 1 N 2 {NN p + N 1 N 2 ) 



(dT/dsY 



N 2 (TV 2 - iy 



(63) 



When Ni — N 2 — 1 we reproduce Eqs. (H0J) and (HTJ. 

In this case, ((dT/dq) 2 ) does not diverges for j p = N p = 




FIG. 3: (dT/de) 2 as a function of 7 P = T v Np for an asym- 
metric cavity connected to two leads with Ni = N2 = 2 open 
channels. The errorbars indicates the result of numerical sim- 
ulations for N p = 200 and T p = 0.025, 0.05, 0.075, 0.1, 
0.125, and 0.15, that give j p = 5, 10, 15, 20, 25, 30. The 
continuous line is the analytical formula given by Eq. 11471 for 
a) (3=1, and Eq. (SSJ for b) (3 = 2. 



0, in contrast with the (3=1 case. Also, this agree with 
Ref. 

Although 7 p takes only integer values, a simple extrap- 
olation to non integer values works qualitatively well as 
can be seen in Fig. 0b), where we have compared Eq. 
(|63|l for Ni = N 2 = 2, with the results from numerical 
simulations 0] for T p = 0.025, 0.05, 0.075, 0.1, 0.125, 
and 0.15 with N p = 200. 

Similarly, we substitute Eqs. (|58|l and (55} into Eq. 
Eq. (|77|) for (3 = 2 to obtain the fluctuations of the 
derivative of R; we have 



(dR/dx)- 



(dR/de) 2 ) = 



ANUdR/de 

2N?(N - iVi) 2 
N *(N 2 - l) 2 



(64) 
(65) 



In contrast with the (3 = 1 case, ((dR/dq) 2 ) does diverges 
at all for (3 = 2. 



IV. FLUCTUATIONS OF dT/dq AND dR/dq 
(q = e, x) FOR SYMMETRIC CAVITIES 

Because of the left-right symmetry of the cavity it is 
sufficient to consider dT/dq, the results for dR/dq are 
equivalent. Also, as happens in asymmetric cavities, 
((dT/dx) 2 ) is always 4N times ((dT/de) 2 ) [see Eqs. O 
and Q62[)]. Then, we will concentrate on the variance of 
the energy derivative of T. 

For LR-symmetric cavities we define o' ab as the 
channel-channel transmission probability, i.e. the square 
modulus of each element t' ab of the transmission matrix 
t' of Eq. (|16fl . It can be written as 



[(o-i)ab + (o- 2 )ab - 2Re/ afc ] , 



(66) 
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where the prime on the left hand side indicates that it is 
defined for LR-symmetric cavities, while ai, cr 2 are de- 
fined by Eq. (|21[) and correspond to Si and S 2 matrices; 
fab is an interference term given by 



fab = {Sl)ab (S 2 ) a fc ■ 

The energy derivative of T is given by 



,6=1 



dT/de = ]T da' ab /ds 
and its fluctuation by 



Ni Ni 

OT/d/ef) = E E W$ V 

a, 6—1 a' ,b' — l 



(67) 



(68) 



(69) 



where, analogous to Eq. (|25|) for q — e, we have defined 
the correlation coefficient for the symmetric case as 



n (P) ab - I da 'ab d(J 'a'b> 
u e a'b' — 



de de 

Using Eq. (|Bl)jl we write Eq. (f70j) as 



(70) 



n(/5) ah -I 

a'b' — ~ 



\ a 6 



° E a'b' + n( - r E a'b' 



(71) 



where C s a>b , is given by Eq. 1)41) [1 for (3=1 and Eq. 
J29 for j3 = 2, with AT replaced by TV' = JVi + AT p /2, 
while 



F (/3) aft 

Q ' b ' ^ de de 



dfab df*. 



To arrive at Eq. I|71[) we used the fact that 
Si and S 2 are statistically uncorrelated, equally and 

(0)O-b ^y(8) a ^ 

uniformly distributed such that C 2e a 'V = C^', b ,, 



that we define as C 



e a'b' 



Also, we use the re- 



sults ([d(ai) a b/de][d(a 2 ) a 'b' /de]) = from one side, 
{[d(<J 3 )ab/de\(df a , h '/de)) = (3 = 1, 2) for the other 
side, and finally {(df a b/de)(df a 'b'/de)} = that are eas- 
ily to verify. 



ab 



In order to calculate a , b , we write it explicitly in 



terms of Si, S2 ; it is given by 



r e a'b' — z 



{{Sl)ab{Sl) a 'b' 



d{S* 2 ) ab d(S 2 ) a 'b' 



(Sl) a( 



d(St)a'. 

• de 



Oe 
(Saja' 



de 

d{S* 2 )ab 

de 



.(73) 



The second line of the Eq. (|73|l is zero as was shown in 
Ref.lU 



A. The P = 1 symmetry 



The correlator D 



From the appendix in Ref. |25j, for j3 = 1 we have 



((Sl) a b(Si) a /f,. 

9(S 2 *) ah d(S 2 )^ 



(9c 



+ + (74) 



N'(N' + 1) 



a=l 



(75) 



where we have replaced A" by AT' = N\ + N p /2 and X 
by e. Then, after we substitute Eqs. l|71jl. f75)l and fT5f> 
into Eq. (|73[) we perform the average over the unitary 
matrix W to arrive to the result 



p(i) ab 



(AT' + l) 2 



(76) 



Eq. O witn AT replaced by A/"' gives (ti) = 1/AT' by 
direct integration, together with Eq. (|3"T|) lead us to the 
result 



_ 2(jV' 2 +iV' + 2)( ( 5g'^+^^') 2 , , 
a ' b ' (A^'-2)A^' 2 (Af' + l) 3 ' 



Finally, Eq. gUJ) with AT' instead of AT and Eq. O 

fi) ab 

gives the result for Z?e ^ a , b , [see Eq. (|7Tf)] . 

As for the asymmetric case, several cases are of par- 
ticular interest. A first correlated case is obtained when 
all indices are equal, which gives the variance of the en- 
ergy derivative of the transmission probability between 
two channels symmetrically located, a' aa ; it is 



(72) ((da'aa/def) = 



N'(N' 2 - 1) + (AT' + 3)(N' 2 + N 1 + 2) 

(AT' - 2)N ,2 {N' + 1) 3 {N' + 3) ' 

(78) 

A second correlated case is obtained for a' = a and b' = b 
but a 6, which gives the energy derivative variance of 
the transmission coefficient a' ab between two channels not 
located in a symmetric way; we have 



{do'aJdef) = 



[{N 1 + 1) + (AT' + 2,)]{N' 2 + AT' + 2) 

A{N' ~2)N' 2 {N' + lf(N' + Z) 

(79) 

The last two equations are different because of the reflec- 
tion symmetry of the cavity. At level of the matrices Si 
and S 2 [see Eq. (|1T|) ] . the diagonal elements represent re- 
flection amplitudes, while the off-diagonal ones represent 
transmission amplitudes. In fact, the first term on the 
right hand side of Eqs. Ij78(l and Q79J1 are equal to Eqs. 
(|49|l and (|41l) (except by a constant factor), respectively, 
when N\ = 1 and N is replaced by N 1 . The second term 
of Eqs. (|78|l and l|79|) comes from interference between 
Si and S 2 [see Eq. (|7TJ) ]. 

In the limit of strong absorption, ({da' a a/de) 2 ) 
and ((do-' ab /de) 2 ) behave as 7I"" 4 . In similar way, 



it is simple to verify that ((da' aa /de)(da' a , b , /de)) 
and ((dcr' ab /de)(da' ab , /de)) behave as 7~ 5 , while 
((da' aa /de)(da' a , a ,/de)), ((da' a Jde)(da' a , v /de)), and 
((da' ab /de)(da' a , b , /de)) go as 7 p 6 . As happens in the 
asymmtric case, the variables da' ab /dq, for a, b = 
1, . . . , Nx, are uncorrelated for strong absorption. They 
enter in the construction of dT/dq [see Eq. (|68H ]. the 
distribution of which is easily obtained when the one for 
da' ab /dq is known [30|. 



jB. Variance ofdT/de 

From Eqs. igUJ) with N replaced by TV', 0, O 
and l|69|) for /3 = 1 we obtain the variance of the energy 
derivative of T, the result is 



(dT/de) 2 ^ = 



Nx(Nx + l) 



2(N' -2)N' 2 (N' + 1) 2 



A + A' + 2 



A' 



1 



(A'- Nx)(N' - Nx + 1) 
A' + 3 



(80) 



The effect of the LR-symmetry is clear. The second term 
of the last equation is similar to Eq. (|49|l with R replaced 
by T. That is because dT/dq for LR-symmetric cavity 
has a similar expression to dR/dq for asymmetric cavity 
as can be seen by comparison of Eq. I|68[) with Eq. (|23H . 
The second term in Eqs. (|80|l comes from the interference 
term of matrices Si, S 2 as explained above [see Eq. J7TJ]. 

For Ni = 1, N' = 1 + N p /2, Eq. (gj reduces to Eq. 
GHJl. In this case ((dT/dq) 2 ) diverges for 7p = N p < 2, 



but remains finite for 7 P 



Nr> > 2. 



((dT/dq) ) diverges only for j p = N p 



When Ai = 2 
0. In both cases 



a complete agreement with the results of Ref. |25| is found. 
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FIG. 4: The same as Fig. [3] but for symmetric cavities. The 
continuous line is the analytical formula given by Eq. I8U1 for 
a) p = 1, and Eq. JUJ for b) /3 = 2. 



Finally, Eqs. JSHJ with A replaced by A', and (JEHJ) gives 
the desired result for D {2)a a , b , [Eq. GU for /3 = 2] . 

In this (3 = 2 symmetry there is not difference in 
the variance of the energy derivative of channel-channel 
transmission coefficient whether the two single channels 
are located symmetrically or not. It is given by 



(da' ab /dey 



1 



A" 



1 



4A' 2 (A' + l) 2 ' 4A' 4 (A' 2 - 1) 



/4/A772 



(84) 



The first term on the right hand side is the same, ex- 
cept by a constant, as Eq. ||BTJ, replacing A by A'. The 
second term comes from interference between Si and S2 
[see Eq. l(7T)l]. For strong absorption, ((da' aa /de) 2 ) be- 
haves as 7~ 3 . Also, as 7 P = N p increases the quantities 
dcr' ab /de for a, b = 1, . . . , Nx, become uncorrelated. 



B. The (3 = 2 symmetry 

1. Correlations of da' ab /dq 

A gain , making an appropriate correspondence from 
Rcf. 25 we have 



((SiUiSZU,) = 5 a a 'S b b /N' 



(81) 



( d(s* 2 ) a b d(s 2 ) a , b> 

\ de 



ds 



^£<(^U-(82) 

a— 1 



We substitute Eqs. JH3 and JT3J into Eq. JZ2| for 

[3 = 2, and perform the average over the unitary matrix 
W, the result is 



(2) «6 _ 2(N' 2 + l)5j8 b b ' 
e a ' b ' N' 4 (N' 2 -1) ' 



(83) 



where we used Eq. (|5fi|l and the result (jx) = 1/A' which 
can be obtained by direct integration from Eq. Q14JI. 



2. Variance ofdT/de 

From Eqs. ^ with A replaced by A', 1^, GB for 
(3 = 2, and 183fl we obtain 



(dT/de) 



N 2 (N'-N!) 2 N 2 (N' 2 + l) 



4A' 2 (A' 2 



I) 5 



4 A' 4 (A' 2 



1) 



(85) 



Again, we note the effect of the LR-symmetry. The first 
term is similar to Eq. I|65l) . The second term in Eq. I|85() 
comes from the interference term of matrices Sx, S 2 - 

For the one channel case (Ai = 1), ((dT/de) 2 ) diverges 
for 7 p = A p = 0, in contrast with the asymmetric case 
for (3 = 2, and in agreement with Ref. I25L It remains 
finite for 7 P > 0. 

In Fig. 2] we compare the analytical results (|80|l and 
(|85p. obtained with T p = 1, with the results from numer- 
ical simulations for T p < 1; we observe a good qualitative 
agreement. 
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C. TRI broken by a magnetic field 

When TRI is broken by a magnetic field, the problem 
of a LR-symmetric cavity is reduced to the problem of 
asymmetric cavity with 0=1 symmetry but the roles of 
T and R interchanged, such that the parametric deriva- 
tive of T is given by Eq. Ij23|l . All the elements da a b/dq, 
for a, b = 1, . . . , Nx, are uncorrelated in the strong ab- 
sorption limit. 

In this case, for instance, the variance of dT/de is given 
by 

/mr/8 ) 2 \ = ^(m + lKN-NjjN-N. + l) 
V 1 ' I (N-2)N 2 (N+l) 2 (N + 3) 

(86) 

For a cavity connected to two leads each one support- 
ing one open channel, ((dT/de) 2 } diverges for j p — N p = 
0, also in contrast with the (3 = 2 case for asymmetric 
cavities. 



V. SUMMARY AND CONCLUSIONS 

The purpose of the present paper was the study of the 
statistical fluctuations of the derivative of the transmis- 
sion T and reflection R coefficients, with respect to the 
incident energy E and an external parameter X (shape 
of the cavity for instance), for ballistic chaotic cavities 
with absorption. 

Our analytical results were obtained assuming N p 
equivalent absorbing channels that are perfectly coupled 
to the cavity (T p = 1). This restrict our calculations to be 
valid in the strong absorption limit, and the absorption 
strength takes only integer values (-f p = N p ). However, 
the results presented here are also valid for no absorp- 
tion, which means 7 P = N p = 0; they are in complete 
agreement with those obtained from known distributions 
of the parametric derivatives of T and R existing in the 
literature. Also, we have shown, by comparison with nu- 
merical simulations, that a simple extrapolation to non 
integer values of -f p is qualitatively correct. 

We considered both asymmetric and left-right (LR) 
symmetric cavities connected to two waveguides: ATj 
channels on the left and N 2 channels on the right; both 
symmetries, the presence and absence of time-reversal 
invariance (TRI), were analyzed. For all cases, the fluc- 
tuations of the energy derivative are smaller than those 
with respect to parametric. We found that ((dT/dx) 2 ) = 
AM '((dT/de) 2 }, where e = 2nE/A and x = X/X c with 
A the mean level spacing and X c a typical scale for X. 
M = N for asymmetric cavities, with N — N% + N2 + N p , 
while M — N/2 for the symmetric case (Ni = N 2 ). 

The correlation coefficient for the parametric deriva- 
tive of the channel-channel transmission probability a a b, 
d(Jab/dq (q = e, x), was calculated. It was shown that 
in the strong absorption limit the different quantities 
d(J a b/dq for become uncorrelated variables. They enter in 
the construction of dT/dq. This is a relevant simplifica- 



tion when the distribution P(dT/dq) is desired assuming 
the one for da a b/dq is known. That is the case of Ref. 
I30I where numerical simulations show evidence of an ex- 
ponential decay for P(d<r a b/de). The decay constant A 
can be obtained directly from ((dT a b/de) 2 ) = 2/ A 2 . A 
similar behaviorfor da a b/dx is expected. This is in con- 
trast with the case of zero absorption where a long tail 
distribution is obtained for the parametric conductance 
velocity 0H|. 

In the case of an asymmetric cavity connected to two 
leads each one with one open channel (N\ = N2 = 1), 
at zero absorption, we find that ((dT/dq) 2 ) (q = E, X) 
is finite when no TRI is present, but is infinite in the 
presence of TRI, in agreement with Ref. |3 where a long 
tails distribution for dT/dq was obtained. The diver- 
gence in the second moment is suppressed by absorption 
and we expect that the long tails become exponential at 
sufficiently large 7 P as mentioned above. This case also 
corresponds to one of an asymmetric cavity with onc- 
lead-one-channel (N\ = 1, JVjj = 0) with one channel of 
absorption perfectly coupled to the cavity, i.e 7 P = 1. In 
this case, ((dR/dq) 2 ) is infinite (finite) in the presence 
(absence) of TRI. ((dR/dq) 2 ) = at zero absorption, as 
should be, and it is infinite for < 7 P < 1. The diver- 
gence disappear for j p > 1. 

For a left-right (LR)-symmetric cavity connected to 
two waveguides with one open channel each one (Ni — 
N 2 = 1), ((dT/dq) 2 ) is divergent for < 7 P < 2, and 
remains finite for "f p > 2 in the presence of TRI. In the 
absence of TRI, the results are different in the presence 
or absence of an applied magnetic field. However, in both 
cases ((dT/dq) 2 ) diverges at j p = 0, in contrast to the 
asymmetric case, and in agreement with Ref. l25l a long 
tails distribution for dT/dq was found at zero absorption 
for presence and absence of TRI. ((dT/dq) 2 ) is finite for 
7 P > 0. We also expect that the long tails will be sup- 
pressed at sufficiently strong absorption 30]. 

The results obtained in this paper help to understand 
some results presented in Ref. |30j about the energy 
derivative of the transmission coefficient, and can serve 
as a motivation to extend that analysis to study the dis- 
tribution of the transmission derivative with respect to 
shape deformations, as well as to motivate the analysis 
of the distribution of the parametric derivative of the re- 
flection coefficient. 
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APPENDIX A: THE COEFFICIENTS M(a,f3,j,6) 

Applying the result (6.3) of Ref. to our case, we 
can write Eq. (|30l) as 

J(a, (3, 7, S) = Ait! + Bu 2 + Cu 3 + Du 4 + Eu 5 , (Al) 
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where 



.4 



B 



iV 4 - 87V 2 + 6 



C = 



D 



N 2 (N 2 - 


1){N 2 -A)(N 2 - 


9) 




N(N 2 - 4) 




N 2 (N 2 


- l)(N 2 -4){N 2 


-9) 




2N 2 - 3 




N 2 (N 2 - 


l)(iV 2 -4)(iV 2 - 


9) 




N 2 + 6 




N 2 (N 2 - 


1)(N 2 -4)(N 2 - 


9) 




5N 





(A2) 



and 



E ' N 2 (N 2 - 1)(N 2 - i)(N 2 - 9)' 

+ ^(8^8*8^,) + ^(8^,8*) 

+ a 5 (8^8*8^8*,) + ae(8^8*8*6^,) 

+ a r (8*8f8^8*) + ^(8*8^4,) 

+ 00(8*6$ 8«,8*,) + a 10 (5*5$6Z,5%,) 

+ a n (8*S*S^,)+a 12 (8*8l'8^,) : 



(A3) 



with 



«3 



1 4 
(1 



6 K K 



€ a 2 = (1 

-SfsM a, = (l + 8^,)8 b 

«5 = (S a b ' + #<#) W a 6 = (S b b ' + 6S'6«)6 b a 5Z 

a? = (5i + «) « a 8 = (8* 4- 8°' 8*)8%8 b , 

(1+8*8^ a 10 = (l + 



(A4) 



Of) 



«n = (1 + S b a 5 a b )8 a a S b ' a 12 = (1 + 8% 5 b b )S b a 8% 



b'\Xb' i 



The coefficients Uj, for j = 2, ... ,5, are obtained from 
u\ through appropriate place permutations of the upper 
indices (a,0,c',c') of the coefficient M of Eq. (|3*U)l . u 2 
is obtained by the sum of the place permutations (12), 
(13), (14), (23), (24), (34), while u 3 by the sum of the per- 
mutations (123), (132), (124), (142), (134), (143), (234), 
(243); u 4 by permutations (12) (34), (13) (24), (14) (23), 
and finally u$ by the place permutations (1234), (1243), 
(1324), (1342), (1423), (1432). The results for u 2 , u 3 , u 4 , 
its are of the same form as Eq. i|A3|l but with a k replaced 
by coefficients that we call b k , Ck, d k , e k , respectively; 
they depend on sums of a k 's. We will see below that not 
all them contribute to Eq. (|29(l : then, we show only the 
coefficients indexed by k — 11, 12 that are important to 
that equation: 



bn 
b 12 

Cll 
Cl2 

du 
di 2 



a 3 - 
a 2 - 
a 2 - 
en 
oi - 
oi - 



a 5 
a 6 
«3 

a 4 
a 4 



a 7 4 
as 4 

an 
an 



a 9 + an + a 12 

ttio 4- an 4- ai 2 

a 6 + a 7 4- a s 4- a 9 4- aio 



(A5) 



en = ai + a 2 4- a 4 4- a 6 + a 7 4- aio 
ei2 = ai + a 3 4- a 4 4- a 5 + a 8 4- a 9 . 



For instance, the result for J(a, /3, 7, <5) can be written 



as 



J(a,P,j,8) 



mi (8«8$8*8*)^ 
m 3 (5°S*5*4) - 
m,(8^8*8 a a ,8*,) - 
m 7 (8* 8^8^,8*,) - 
m g (8*8$8^8*)- 
m n (8*8*SZ4,) 



-m 2 (8^8*8^8*,) 
-m 4 (S^8f8*S*,) 
-m 6 (8^8*8*8^) 
-m 8 (8*8%8i,sP,) 
-m 10 (S*8$Si,Sp) 
+ m 12 (6*5l'5?,6%,), 
(A6) 



where 



m k = Aa k + Bb k + Cc k + Dd k + Ee k , k = 1, . . . , 12. 

(A7) 

From Eq. (| A6|> we construct the coefficients J(a, (3, c, c) 
and J(c, c, a, (3) , take the difference of them and sum with 
respect to c, d . The result is given by Eq. H31|l. where 



m = mn 4- mi2 

n 2 — m% — m 3 — mg + mio — Nmn 
7i 3 = m 2 — m 3 — mg + mio + Nmi 2 . 

From Eqs. (|A9|) and (|A10|) we see that 

"3 - "-2 = ATrai 



(A8) 
(A9) 
(AIO) 



(All) 



Eqs. and (J3TJ leads to Eq. (J22Jl, the result being 
dependent on n\ , and n 2 , 77,3 through the difference n 3 — 
719 = Nni. From Eqs. (|A5|I . (|A7|I and (|A8|I . «i is given 
by 

ni = (A + 2 J B + Z))(aii4-ai2)4-2( J D + J B)(ai+a 4 ) 
+ (S4-2C + S)(a 2 4-a 3 4-a 5 
4- a 6 + a 7 4- a s 4- a 9 + aio). (A12) 



Finally, Eqs. I|A2|I and l|A4jl give 
1 



Tlx 



N 2 (N 2 - l)(N + 2)(N + 'i) 
+ (N + l)(N + 2)(5 a a 'S b b +5*5°') 

- (N + l) 



{2(1 + ^)(1 + ^,) 



°b 



+26 b J* (d b '6 a a , + St'sp) + 2 (6*81,8%: 

+ 5 b a 6 a b '6 a a , + S b a 5 b 'S a b , + 8*81,8$)] } (A13) 
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